Distinguishing quantum and classical transport through nanostructures 
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We consider the question of how to distinguish quantum from classical transport through nanos- 
tructures. To address this issue we have derived two inequalities for temporal correlations in nonequi- 
librium transport in nanostructures weakly coupled to leads. The first inequality concerns local 
charge measurements and is of general validity; the second concerns the current flow through the 
device and is relevant for double quantum dots. Violation of either of these inequalities indicates 
that physics beyond that of a classical Markovian model is occurring in the nanostructure. 



Quantum coherence of electrons is the essential ingre- 
dient behind many interesting phenomena in nanostruc- 
tures (e.g., 0, Q). Considerable progress has recently 
been made in the investigation of coherent effects in 
nanostructures with both charge and transport measure- 
ments, (e.g., 0-0] )■ Typically Rabi oscillations in the 
current are taken as a distinctive signature of quantum 
coherence. However, since even classical autonomous rate 
equations can admit oscillatory solutions (e.g., d), oscil- 
lations by themselves cannot be considered as a definitive 
proof of the existence of quantum coherent dynamics. 

In this paper we formulate a set of inequalities that 
would allow an experimentalist to exclude the possi- 
bility of a classical description of transport through a 
nanostructure. The inspiration for this comes from the 
Leggett-Garg inequality Q, which has been described as 
a single-system temporal version of the famous Bell in- 
equality, also a topic of interest in nanostructures at the 
present (e.g., (H). The Leggett-Garg inequality @ can 
be summarized as follows. Given an observable Q(t), 
which is bound above and below by \Q(t)\ < 1, the as- 
sumption of: (Al) macroscopic realism and (A2) non- 
invasive measurement implies the inequality, 

(Q(h)Q) + (Q(h+t 2 )Q(ti)) - (Q(h+t 2 )Q) < l, (l) 

where Q = Q(t = 0). The question of (Al) 'realism' Q 
can be phrased as: before we perform the measurement 
Q on the system [llj, does it have a well defined value? 
A classical system does, but a quantum system does not. 

In the context of nanostructures weakly coupled to 
contacts, such that a generalized master equation de- 
scription (e.g., (EE El) is appropriate, we derive and 
study two inequalities. The first concerns correlations 
between local charge measurements performed, e.g., by 
a quantum point contact (QPC) (e.g., Q)- We formu- 
late this inequality in quite general terms, applicable to 
a range of nanostructures. The charge measurements 
we consider here are related, in spirit, to recent work 
(e.g., on violations of the Leggett-Garg inequality, 



CD W77777777, 




V////////A 



using continuous weak measurements on closed systems. 
However, in contrast to their work, here we are consider- 



FIG. 1. (Color online) (a) shows a generic single-charge 
nanostructure. The single-charge can occupy one of N in- 
ternal states. Local charge occupation of one or more sites is 
measured using a charge detector (CD) Eq. ([2}, e.g. a quan- 
tum point contact, (b) shows a schematic of a double quan- 
tum dot, the main example we discuss here. F L / R describe 
the left /right tunnelling rates, A is the coherent tunnelling 
amplitude between the left and right dots, and e is the energy 
difference between the left and right occupation states. Local 
charge measurements are done on the right (R) state for the 
charge inequality of Eq. ([2]), or transport is measured into the 
collector reservoir for the current inequality of Eq. (|4|). 



ing a very different situation: ensemble averages of strong 
(i.e., projective) non-continuous measurements on open 
transport systems. Moreover, our second inequality ex- 
plicitly focuses on DQDs, providing an inequality for the 
correlation functions of the current flowing through this 
widely-studied nanostructure. This second inequality is 
of particular relevance to DQD experiments along the 
lines of those of Ref . 0, Hj , where we predict that viola- 
tions of both inequalities should occur. 

Systems. — We begin by outlining the class of sys- 
tems studied in this work. We consider nanostructures 
(Fig. Q]) weakly coupled to leads such that transport pro- 
ceeds via sequential tunnelling, and we assume a large 
bias such that higher-order tunnelling, level-broadening, 
and non-Markovian effects can be neglected [HI]. We 
assume strong Coulomb blockade such that the system 
admits at most one excess electron. In these limits the 
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master equation formalism we apply here, while simple, 
has been shown to be very accurate by a variety of experi- 
ments 0-0| • In general, non-Markovian effects might lead 
to a violation of these inequalities, so care must be taken 
to verify one is in these limits. Our system comprises 
of (N + 1) states: the "empty" state, |0), with no ex- 
cess electron, and states \n) with a single excess electron 
in state n = 1 . . . N. The dynamics is described by the 
generalized master equation p(t) = W[p(t)], where p(t) 
is the density matrix and the superoperator W is the Li- 
ouvillian. Within a quantum-mechanical description, the 
density matrix p(t) contains coherences, and the Liouvil- 
lian has a Lindblad form W = Wq+T,, where the coherent 



evolution of the system is given by Wq 
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H the internal system Hamiltonian, E is the self-energy 
induced by the contact with reservoirs (both electronic 
and otherwise), and throughout we set e = h = 1. 

The analogous classical description is a rate equation 
for the probabilities P n (t) of finding the system in state n 
at time t. This rate equation can be written in the same 
form as above: p = W cl [/o], but now the density matrix p 
only includes diagonal elements [the probabilities P n (t) = 
Pnn(t)\- This p can be represented as a vector such that 
the Liouvillian W cl is a classical rate matrix with Wf] > 
0;i / j and W^ 1 = — Yljjti^jh O ur general strategy 
is to explore the behavior allowed by this classical rate 
equation and use this to derive our inequalities. 

Charge Inequality. — The first inequality we derive is 
for localized state measurements. Consider a charge de- 
tector which registers the value Q n > when the system 
is in state n (see Fig. [T]). Furthermore, let us designate 
as state N the state for which Q has maximum value: 

Qn Qm&x- 

(i) Classical regime: We assume that for a classical 
system the charge measurement can be performed non- 
invasively. An initial state, described by a set of prob- 
abilities P„(0), is fixed and known (actually one only 
requires knowledge of the relevant expectation values in 
this state) . When non-invasively measuring the charge of 
a classical Markovian system, we posit that the following 
inequality holds 

\L Q (t)\ = \2(Q(t)Q) - (Q(2t)Q)\ < Q max (Q), (2) 

where (Q) — Y^k^i^Qk is the expectation value of 
Q = Q(t = 0), and (Q(t)Q) is the charge-charge correla- 
tion function. This inequality holds in two regimes: (i) 
stationarity, where it follows from the original Leggett- 
Garg [Eq.(l)] by defining the normalized operator Q — 
2Q/Qmax — 1 [15| and taking the stationary expecta- 
tion value; and (ii) if only a single state contributes 
to the detection process, i.e., Q n — Qm&x^nN, then 
Eq. (|2|) holds for an arbitrary initial state (defined 
by the set of probabilities Pfe), and not just the sta- 
tionary state. This latter can be seen as follows. 
Within both classical and quantum stochastic theory, 
the charge-charge correlation function can be written as 



(Q(t)Q) = Qmax^A'Ar(OQmaxPv(0), where the "propa- 
gator" r2jvjv(i) is an element of the stochastic matrix giv- 
ing the probability of finding the system in local charge 
state N a time t after it is in state N. The quantity 
Lq can thus be written as (for this single state mea- 
surement) L Q (t) = Q 2 max p N (o)[2n NN (t)-n NN (2t)}. 

If the behavior is classical and Markovian, then the 
Chapman-Kolgomorov equation for classical rate equa- 
tions applies 13(| , and we can write the propagator 



with argument 2t as a decomposition over intermediate 
states fijvjv(2i) = Sfe ^Nk(t)^kN(t) to obtain Lq(£) = 
Q"La. x p N(0)[^NN{2 - Qnn) - Ysk^N^Nk^kN}- Here- 
after, we suppress the time argument, Q — f2(t). £q(£) 
is then maximized by choosing the propagators such that 
the system always ends up in state N, i.e. £Inn = 1, 
which gives: max{ L Q (t)} = Q, 2 nax Pv(0) = Qmax(<3)- 
The lower bound is: min{L<g(t)} = — Q ma .x(Q)- For this 
single state measurement the inequality holds indepen- 
dent of initial state, as the dynamics are sufficiently con- 
strained by the Chapman-Kolgomorov equation alone. 
We first illustrate the charge inequality violation with 
an example, before continuing to derive the inequality 
for current measurements. 

(ii) Quantum regime: The transport DQD consists of 
a dot L, connected to the emitter, and dot R, connected 
to the collector (see Fig. [T]). Assuming weak coupling, 
large bias, and Coulomb Blockade, the basis of electron 
states is {|0), \L), \R)}. Its Hamiltonian becomes 

H = e(\L)(L\ \R)(R\) + A(\L)(R\ + \R)(L\), (3) 

with e the level splitting, and A the coherent tunnelling 
amplitude betweens the dots, and with self-energy 

- E 

a=L,R 

where sl = |0)(L|, sr = |i?){0|, and Tl and Tr are the 
left /right tunnelling rates (throughout we set e = H = 1). 
The influence of phonons can also be included in £ in 
the standard way [l|, [l6j]. The corresponding classical 
Liouvillian is a 3 x 3 matrix with elements Wr^; a, (3 — 
0,L,R. For illustrative purposes, we consider a charge 
measurement in which the detector only registers when 
there is an electron in the righthand QD: Q = \R){R\ 
for which Q max = 1. The correlation functions are then 
calculated from (Q(t)Q) = Tr [Qe wt Qp ] , with p the 
stationary density matrix of the system. 

In Fig. [2]we plot \LQ(t)\/(Q max (Q)) as a function of 
time for a DQD. The behavior is oscillatory, but also 
damped due to coupling to both the collector and the 
phonon bath. The shaded region (> 1) indicates where 
|Lq(£)| violates the inequality of <j2j> . The most promi- 
nent violation occurs at the maximum closest to t = 0. 
For these parameters then, no classical Markov descrip- 
tion of the system is possible and here, it is quantum 
oscillations between L and R states that are responsi- 
ble for the violation. As we discuss later, the degree of 
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FIG. 2. (Color online) The charge inequality 

\LQ(t)\/(Q max (Q}) Eq. © for QPC charge measure- 
ments on the R state of a double quantum dot, as a function 
of dimensionless time At. The parameters are shown on the 
figure, and we have set Q max = 1- We chose a dimensionless 
coupling constant g — 0.05 for the bulk phonons 1]. The 
solid line represents the case of no phonons, and the dashed 
line includes a phonon bath at temperature T = WhA/ks, 
to illustrate how an invasive environment masks the viola- 
tion. The colored region marks the area of violation of the 
inequality in Eq. @. 

violation can be increased by decreasing T^, which per- 
mits the electron to spend a longer time in the DQD. 
In the limit T R — > and with Tl — > oo, such that 
the empty state may be eliminated, we find the ana- 
lytic form 7 Q (t)/(Q max (Q)) = [cos(2At) + sin 2 (2At)] , 
for e = 0. The coherent tunnelling A defines the time 
when the violation is maximum, t max = 7r/6A, such that 
LQ(t max )/(Q max (Q)) = §. The time t max is in agreement 
with that observed for the Leggett-Garg inequality [9( 
for a single free qubit with level coupling A/2. The ef- 
fects of a phonon bath are also apparent in Fig. [5J where 
we have used reasonable bath parameters Although 
the oscillations of 7g(t) are damped, the first and most 
significant maximum remains. 

Current Inequality. — Our second inequality concerns 
the current I(t) flowing through the transport DQD: 

|L / (t)| = |2(/(t)7)-(7(2t)7)|<r ii (7) ) (4) 

where T R is the coupling to the collector, I = I(t = 0) 
and (7) is the average current of the initial state. Al- 
though this second inequality resembles the first one, 
in Eq. @ (T R is the maximum instantaneous collec- 
tor current), its derivation and significance are some- 
what different. This is because, in the master equa- 
tion approach, the current operator translates into a 
"jump" super-operator and Eq. (UJ) thus represents an 
inequality concerning transitions in the system, and 
not static properties such as the charge under the 
noninvasive measurements of Eq. ([2"]l. For the DQD 
model in the infinite bias limit, the current super- 



FIG. 3. (Color online) The current inequality 

|7j (£)l/(rj?{7i{)) Eq. for current measurements on a dou- 
ble quantum dot, as a function of dimensionless time At. The 
solid line represents the case of no phonons, and the dashed 
line includes a phonon bath at T — lOfiA/fcs. The colored 
region marks the area of violation of the inequality in Eq. Q . 
The violation of this inequality is more pronounced than that 
of the charge inequality. The effect, however, is more sensitive 
to the phonon bath. 

operator acts as J[p\ = T R \0)(R\p\R)(0\, such that 
the average current is (7) = Tr { J p} and the cor- 
relation function of interest is obtained as (I(t)I) = 
Tv{je wt Jp a }, where again the stationary distribution 
is chosen as the initial state. In these terms, Eq. (|4| 
can be written as 7/(t) = Tr {j (2e m - e 2Wt ) Jp Q } . 
In the classical description of the DQD, J is the 3 x 
3 matrix with elements J a j3 = T R S a fi5i3 tR . Thus 
using Chapman-Kolgomorov again, we have 7/(t) = 
T 2 R P R (0) (fl R0 (2 - n 00 - n RR ) - n RL n L0 ) . For a gen- 
eral Markov stochastic matrix, f2, the maximum of 7/ is 
2T R P R (0). However, the rate equation form Q(t) = e wt 
furnishes us with a further constraint. Maximizing 7j(t) 
with respect to time, from 7/ = and SI — Wf2, we find 
that the maximum of 7/ occurs when fioo + Sl RR = 1 and 
fl R0 = l, giving max{7 / } = T 2 R P R {0) = T R (I). This re- 
sult relies on the geometry of the DQD, and in particular 
the form of the jump operator and the absence of direct 
tunnelling from emitter to dot R, i.e. W R o = 0. 

Figure [3] illustrates the violation of the current in- 
equality Eq. Q for the DQD. As with the charge mea- 
surement, the quantity 7/(t) is oscillatory and violates 
the respective inequality with the strongest violation oc- 
curring at the first maximum, which is here at a time 
i max = 7r/(2A). The degree to which this current in- 
equality is violated is of greater magnitude than that for 
the charge measurement. Again, in the limit — > oo, 
one can eliminate the empty state and find an ana- 
lytical form. In addition, the T R — > limit gives 
Li(t)/(T R (I)) = -2sin 2 (At)cos(2At). Thus the viola- 
tion has a maximum of 7/(t max ) / '(T ' R (I)) = 2. 

Under the above assumptions, the three-state classical 
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FIG. 4. (Color online) 

max{Z/Q(i)/Q max (<2)} and (b) 



Contour plot of (a) 
max{Li(t)/T R (I R )} as a 



function of the tunnelling amplitude A/F_r and relative 
rates F_l/F.r. As discussed in the text, violation requires 
A > 0.3r B in (a) and, in (b), A > 1.5F R , T L /T R > 1. The 
phonon bath is excluded. 



DQD Liouvillian cannot produce a violation. However, if 
these assumptions are relaxed (e.g., allowing Wro ^ 0), 
a small violation (on the order of 0.003% Tr(I)) of the 
inequality can be observed in extreme parameter regimes. 
This is in contrast to the Lq inequality, where no further 
constraints are required of the Liouvillian. This differ- 
ence reflects the fact that here the current measurement 
is essentially a destructive measurement of the state of 
the nanostructure. An infinitesimal time interval after 
a positive current measurement is obtained, the electron 
has left the system, leaving it in the 'empty' state. This 
behavior is implicit in the jump super-operator form of 
the current measurement. 

Figure U shows how the maximum degree of violation of 
both inequalities depends on the parameters of the DQD 
with no phonons. A violation of the current inequal- 
ity Eq. gj) requires A > I.STr, Tl S> Fr, and small 
detuning (e < A). The violation of the charge inequal- 
ity Eq. @ is more resilient, and always occurs unless 
there is strong over-damping from the reservoir 3> A. 

Finally, we note that in practice one needs to measure 
the correlation functions in Eq. ([2]) or Eq. (|4|) on very 
short time scales (e.g., 0). Alternatively, one can ob- 
tain either correlator from the inverse Fourier transform 
of the appropriate noise power function. In the trans- 
port case (Eq. Q), one must consider contributions from 
both particle- (as considered here) and displacement- 
currents. In principle, one can either choose appropriate 
gate/junction capacitances to neglect the displacement 
current contribution, or include them in the definition of 
Eq. (Q}, and its subsequent maximization. 

Conclusions. — In summary, we have derived two 
inequalities for non-equilibrium transport in nanostruc- 
tures: one concerning local charge measurements and the 
other for current flow through the device. The first is 
of general validity; the second of relevance to the usual 
DQD geometry found in numerous experiments. Viola- 



tion of either of these inequalities indicates that physics 
beyond that of a classical Markovian model is occurring 
in the nanostructure. This may be taken as evidence for 
quantum oscillations of the electron within the device; 
or it may indicate a non-Markovian interaction with pre- 
viously unappreciated degrees of freedom. Finally, we 
point out that these ideas can be expanded in a number 
of different directions: to other types of measurements; 
to inequalities with different time-dependencies, as in the 
original Leggett-Garg work Q; and to different physical 
situations for which master equations are appropriate, 
such as atom-field interactions in quantum optics. This 
work can also be applied to networks of quantum dots, 
Cooper pair boxes, and molecules. 
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The rates for the electron-phonon interaction are, 
7± = [-e(T/A 2 )Coth[^A/2] ± T/A]7r/2J[A[ 
(T 2 /A 2 ) Coth[/3A/2]27rJ[A], where A 
J[A]:=(?Ae' _A,/nc ', g is the dimensionless coupling 
strength, /3 = (k(,T) _1 and Sl c is an ultraviolet cutoff. 
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